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ABSTRACT 

(N 

Aims. Chemical evolution models are useful for understanding the formation and evolution of stars and galaxies. Model predictions 
will be more robust as more observational constraints are used. We present chemical evolution models for the dwarf irregular galaxy 
NGC 6822 using chemical abundances of old and young Planetary Nebulae (PNe) and H n regions as observational constraints. Two 
sets of chemical abundances, one derived from collisionally excited lines (CELs) and one, from recombination lines (RLs), are used. 
We try to use our models as a tool to discriminate between both procedures for abundance determinations. 

Methods. In our chemical evolution code, the chemical contribution of low and intermediate mass stars is time delayed, while 
for the massive stars the chemical contribution follows the instantaneous recycling approximation. Our models have two main free 
parameters: the mass-loss rate of a well-mixed outflow and the upper mass limit, M up , of the initial mass function (IMF). To reproduce 
the gaseous mass and the present-day O/H value we need to vary the outflow rate and the M up value. 

Results. We calculate two models with different M up values that reproduce adequately the constraints. The abundances of old PNe 
are in agreement with our models and support the star formation history derived independently from photometric data. In addition, by 
assuming a fraction of binaries producing SNIa of 1%, the models fit the Fe/H abundance ratio as derived from A supergiants. The 
first model (M4C), that assumes M up = 40 M Q , fits, within errors smaller than 2 cr, the O/H, Ne/H, S/H, Ar/H and Cl/H abundances 
obtained from CELs, for old and young PNe and H n regions. The second model (MIR), that adopts M up = 80 M Q , reproduces, within 
2 cr errors, the O/H, C/H, Ne/H and S/H abundances adopted from RLs. Model MIR does not provide a good fit to the Cl/H and Ar/H 
ratios, because the SN yields of those elements for m > 40 M Q are not adequate and need to be improved (two sets of yields were 
tried). From these results we are not able to conclude which set of abundances (the one from CELs or the one from RLS) represents 
better the real abundances in the ISM. We discuss the predicted AY/ AO values, finding that the value from model MIR agrees better 
with data for other galaxies from the literature than the value from model M4C. 

Key words, galaxies: individual: NGC 6822 (DDO 209) - ISM: planetary nebulae: general - ISM: Hn regions -galaxies: chemical 
■ evolution 
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Q\ ' 1 . Introduction Timmes et al. 1995; Romano et al. 2010; Kobayashi et al. 201 1), 

*n the authors were not able to test Ne, CI, and Ar yields, due to the 

Planetary Nebulae (PNe) constitute one of the most valuable j ac j c Q f s t e ji ar abundances 
f- } chemical tracers of the past abundances in the interstellar 

^ medium (ISM). Their chemical compositions allow us to deter- The chemical evolution equations (e.g., Tinsley 1974) take 

_ i mine the abundances of some chemical elements present in the into account many physical parameters: galactic infalls, galactic 

^ ■ ISM when their progenitor stars were born. PNe are produced outflows, the initial mass function, the star formation rate, and a 

by stars with initial masses from ~ 1 M to ~ 8 M Q and also set of stellar yields for different masses. Therefore these equa- 

with a large age spread (from 0.1 to 9 Gyr, Allen et al. 1998). tions are complex and have to be solved using numerical meth- 

Therefore PN characteristics are important as observational con- ods. However they can be simplified assuming the instantaneous 

straints in chemical evolution models, allowing us to improve recycling approximation (IRA, Talbot & Arnett 1971). For this 

the inferred chemical history (Hernandez-Martinez et al. 2009, approximation the lifetimes of all stars more massive than 1 M Q 

hereafter HPCG09; Richer & McCall 2007; Buzzoni et al. 2006; are negligible compared with the age of the galaxies. This ap- 

Maciel et al. 2006). Despite the fact that PNe show some bright proximation allows us to solve the chemical evolution equations 

emission lines, deep observations are needed to determine their analytically. Despite its simplicity, IRA is a good first approxi- 

physical conditions and accurate chemical abundances, which mation for elements produced mainly by massive stars (MS), but 

are based on much fainter lines. not for the elements produced partially by low and intermediate 

In addition, gaseous nebulae are an important key in the mass stars (LIMS). 

chemical abundances determination of noble gases (e.g., Ne and There are intermediate methods to calculate chemical evo- 

Ar) and other elements like CI and, therefore, in the test of stel- i ut i on mo dels, which consist in analytical approximations that 

lar yields of these elements. The determination of this type of consider the delays in chemical enrichment produced by LIMS. 

elements in stars is not so reliable and, in previous papers (e. g., Several authors have presented their own analytical approxima- 

tions (e.g., Serrano & Peimbert 1983; Pagel 1989; Franco & 

Send offprint requests to: L. Hernandez-Martinez Carigi 2008). They propose some time-delay prescription for 
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the chemical enrichment produced by LIMS. These time-delay 
terms make the LIMS to bring out to the ISM the processed nu- 
clear material at a single time after their formation, while the 
contribution due to MS is instantaneous, like in the IRA approx- 
imation. 

In this paper we calculate chemical evolution models for 
the dwarf irregular galaxy NGC 6822 following the method 
used by Franco & Carigi (2008). This method was modified by 
Hernandez-Martinez (2010) to include numerically infalls, out- 
flows, and star formation rates. Also, in this new code we have 
increased the number of chemical elements considered, from 5 
to 27. 

NGC 6822, a galaxy of the Local Group, is located at 460 
kpc from the Milky Way (Gieren et al. 2006). It presents a recent 
increase in the star formation rate as shown by its bright H n re- 
gions. These features make it easy to determine the present-day 
chemical abundances of the ISM. Thus, it is suitable for chemi- 
cal evolution modeling. 

Carigi et al. (2006, hereafter CCP06) performed chemical 
and photometric evolution models for NGC 6822. Based on a 
cosmological approach they obtained the gas infall rate adequate 
to form the galaxy and, based on the photometric properties, they 
derived a robust star formation history. Their chemical evolution 
models were built to reproduce the present component of the 
ISM, as given by the chemical abundances of the H n region HV, 
determined from recombination lines (RLs). 

HPCG09 determined abundances from collisionally excited 
lines (CELs) for 1 1 PNe and one H n region. From these data, 
they confirmed the chemical homogeneity of the present com- 
ponent in the ISM and found the presence of two populations 
of PNe. Based on these results they built a preliminary chem- 
ical evolution model to reproduce the chemical behavior of 
NGC 6822. One of the aims of this work is to compute detailed 
chemical evolution models by using the abundances of the ISM 
past component (represented by PNe) as an additional restriction 
to the one posed by H n regions to the ISM current component. 

Moreover, as it is known that abundances obtained from RLs 
are larger that those from CELs, by factors of about 2, a prob- 
lem known as "the abundance discrepancy" (see §5), we also 
compute models to reproduce the chemical abundances as de- 
rived from RLs. Then, by means of these models we try to dis- 
criminate between the values of abundances obtained from both 
methods. 

Also, chemical evolution models can be used to constrain the 
stellar yields, when the model results are compared with accurate 
observations. In this work, we try to fit the observed CI and Ar 
abundances, and derived SN yields for massive stars. 

Furthermore, we also compute the He to O abundance en- 
richment ratio, AY/ AO (in this expression Y and O are the chem- 
ical abundances of He and O, both by mass), which is an impor- 
tant parameter in chemical evolution discussions. 

This paper is organized as follows: In §2 we describe the 
chemical evolution equations used by assuming a time-delay 
approximation for LIMS. The assumptions for the models are 
presented in §3. In §4 and §5 we present the results from our 
models, with observational restrictions obtained from CELs and 
RLs, respectively. In §6 we discuss the evolution of AY vs. AO. 
Finally, in §7 we discuss our results and present our conclusions. 

2. Chemical evolution equations with a delayed 
contribution from LIMS 

Franco & Carigi (2008) proposed an analytical solution for the 
chemical evolution equations, by considering that the enrich- 



ment produced by LIMS is delayed by the lifetime, t, of a rep- 
resentative mass, while the MS enrich instantaneously the ISM. 
This solution was built only for a star formation rate proportional 
to the gas mass in a closed box model. They obtained results very 
similar to those found from numerical models that consider the 
lifetime of each star. Franco & Carigi studied the evolution of 
five elements: H, He, C, N, O, and the metallicity Z. 

Following the prescription given by Franco & Carigi (2008), 
we developed a numerical code which solves the differential 
equations of the evolution of the gas mass M gas (t) (eq. 1), and the 
evolution of the chemical abundance by mass in the ISM for the 
element i, X,-(f) (eq. 2), respectively (Hernandez-Martinez 2010). 

This numerical code can solve 27 chemical species (H, He, 
C, N, O, F, Ne, Na, Mg, Al, Si, P, S, CI, Ar, K, Ca, Sc, Ti, V, 
Cr, Mn, Fe, Co, Ni, Cu, Zn) for which we have a good collection 
of Z dependent integrated yields for LIMS (yu ms ), MS (y ms ) and 
Type la supernovae, SNIa (y Jm - fl ), for 6 initial stellar metallicities 
(Z = 1.0xl0~ 8 , 1.0xl0~ 5 , 1.0xl0~ 4 , 4.0xl0~ 3 , 8.0x10^, and 
2.0xl0~ 2 ). These yields are described in detail in section 3. 

In the delayed approximation the evolution of the gaseous 
mass is given by: 



dM gas {f) 
dt 



-( 1 - R ms )m + R, im .Mt ~ T m ) + f(t)-w(t), ( 1 ) 



here, R ms and R; imj are the masses returned to the ISM by mas- 
sive and low-intermediate mass stars, respectively. r m is the time 
delay of LIMS for their total gas ejection to the ISM, and i^(t), 
/(t) and w(t) are the star formation, accretion and outflow rates 
as function of time, respectively. 

In addition, the evolution of the gaseous mass in element i 

is: 



djX^MgaAf)) 

dt 



-XfrWf) + Ei(t) + Xj (t)f(t) - X™(t)w(t), (2) 



where, x! (t) and X?(t) are the abundances (by mass) of element 
i in the gas accreted and in the gas expelled from the galaxy, 
respectively. Ej(t) is the gas rate, in the element ;, returned to 
the ISM by MS, LIMS, and SNIa that die at a time f, and it is 
calculated by, 

Ei = [R ms Xi(t) + yi, ms ]if/(t) + [Ri ims Xi(t - Ti) + yij ims ]t//(t - t,)+ 

ji,sma iA(f - T snia ), 

here, t, is the time delay when the group of LIMS enrich the 
ISM with the element i (see Franco & Carigi 2008) and T snia is 
the time delay of SNIa to pollute the gas mass of the galaxy. 

In Appendix A we present, for the six initial stellar metal- 
licities and two different upper mass limits of the initial mass 
function, r m , t ; , R mj , R Ums , y ms , y Ums , and y mia for the elements 
considered in this work. 



3. Assumptions of the Chemical Evolution Models 

In the following we describe in detail the assumptions included 
in our models. 

1) We adopt the mass accretion rate proposed by CCP06, 
which was tailored to NGC 6822 from cosmological models of 
galaxy formation. Therefore, the accretion history is a better ap- 
proximation for this galaxy than simple parametric equations. 
CCP06 described two families of cosmological models for the 
increase of baryonic mass: S (small-infall) and L (large-infall). 
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We have chosen the L-models because they adjust better the ob- 
servational constrains presented by CCP06 (the S -models pre- 
dict higher enrichment in O/H and C/H than the L-models). 
In addition, we have parametrized the accretion history of the 
baryonic mass as a function of time, derived by CCP06. Such a 
parametrization, which is used for our chemical evolution mod- 
els, is presented in Fig. Q] (upper panel) and it is given by the 
expression: 

f(t) = 131.4(e-°- 31 '-e-°- 32 ') (10 s M Gyr 1 ). (3) 

The accreted material is assumed to be primordial, that is, 
with a chemical composition given by: X=0.752, 7=0.248, and 
Z = 0.000 as derived from WMAP observations by Dunkley et 
al. (2009) and for metal poor irregular galaxies by Peimbert et 
al. (2007a). 




CO 
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Fig. 1. Accretion (upper panel, see eq. 3) and star formation 
(bottompanel, CCP06) histories adopted in the models. 13.5 Gyr 
corresponds to the present time 



2) We consider the history of the star formation rate (SFR) 
derived for NGC 6822 by CCP06 based on observed color- 
magnitude diagrams and photometric evolution models (Fig. [T| 
bottom panel). Interestingly it is observed that the accretion mass 
history (shown in the upper panel of Fig. [T} and the SFR do not 
follow the same trend and seem to be in contradiction in the lat- 
ter 3 or 4 Gyr, when the accretion diminishes while the SFR in- 
creases abruptly. It should be considered however that the SFR 
should not follow the gas accretion closely as this accretion is 
only one of the mechanisms for increasing the SFR. On large 
scales, the formation of stars is controlled by an interplay be- 
tween self-gravity and supersonic turbulence, operating on the 
gas already accreted. It is not clear how star formation bursts are 
triggered in dwarf galaxies but, for NGC 6822, Gouliermis et al. 
(2010) suggest that turbulence on kpc scales is the major agent 
regulating star formation. Several mechanisms have been pro- 
posed as able to produce turbulence in large scales. In particular, 
for NGC 6822, a possible mechanism, that could have triggered 
at least the increase in the SFR during the last 100-200 Myr, has 
been found in the presence of an apparently independent stellar 
structure called the "northwestern" companion, which is con- 
sidered to be currently interacting with the main body of NGC 
6822 and causing the tidal arms in the southeast of the galaxy 
(de Block & Walter 2000, 2003). If such an interaction occurred 
in the past (which is possible if the northwestern companion is 
a satellite nearby NGC 6822) it could have been the mechanism 
for triggering different starbursts in NGC 6822. 

3) We use the initial mass function (IMF) of Kroupa et al. 
(1993), in the 7.5 - M up M mass range for Z < 10~ 5 and in the 
0.1 - M llp M mass range for Z > 10~ 5 . M up is a free parameter 
adjusted to reproduce the observed O/H abundance ratio. 

4) We assume different stellar yields, dependent on initial 
stellar mass (m) and on initial stellar metallicity Z,-, provided by 
different authors. However they do not cover all the mass range. 
In the case of LIMS, we have a complete library from 1 M to 
6 M and for MS, the mass range covered goes from 9 M to 
80 M . In order to fill the mass gap (from 6M to 9MO), we 
proceeded as follows: for the 6 < m/M Q < 7.5 mass range we 
assigned the stellar yield values given by mass fraction (pi) of 
the 6M star. In the same way, for the 7.5 < m/M < 9 mass 
range we assigned the /?, values of the 9 M star. 

Specifically, 

4. 1) For LIMS we adopt the stellar yields by Karakas 
& Lattanzio (2007), and the integrated yields used by our model 
are given by: 

Yijims = I mpi(m)IMF(m)dm. (4) 

JlM 

The results for ten elements, five initial stellar metallicities, and 
two different M llp values are presented in Tables 4 and 5. 

4.2) For MS, we consider the stellar yields obtained 
from models for the pre-SN stage and the SN stage. The pre-SN 
yields are taken from the work of the Geneva group (Maeder 
1992; Meynet & Maeder 2002; Hirschi et al. 2005; Hirschi 
2007). The SN yields are taken from Woosley & Weaver (1995) 
adopting their models B for the 12-30 M range and their mod- 
els C, for the 35 - 40 M range. We combine the Geneva group 
yields with the Woosley & Weaver yields using the prescription 
proposed by Carigi & Hernandez (2008) where they connect the 
mass of the carbon-oxygen cores (Mco) from the Geneva group 
to Mco from Woosley & Weaver, using the prescription given 
by Portinari et al. (1998). Under these assumptions the He, C, 
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N, and O yields are equal to the pre-SN ones, but for heavier 
elements the yields are similar to the SN ones. For m > 40 M 
the adopted yields for the heavier elements are similar to those 
given for m = 40 M . 

The chemical contribution of MS is in IRA, and the inte- 
grated yields are given by: 

rM„ p 

Yi,ms = I mpi(m)IMF(m)dm. (5) 

Jl.5M a 

The results for ten elements, six initial stellar metallicities, and 
two different M up values are presented in Tables 4 and 5. 

4.3) For SNIa, we consider the stellar yields, P,(snia), 
given in M Q by Nomoto et al. (1997) which are almost indepen- 
dent of the initial stellar metallicity. We assume that a fraction 
Abi„ of the stars with masses between 3 and 15M corresponds 
to binary systems and that every one of these systems becomes a 
SNIa. A/„„ is a free parameter of the model which is found by fit- 
ting the observed Fe/H ratio. It is considered that all SNIa of each 
stellar generation enrich the interstellar medium at T s „, fl =l Gyr 
after the formation of SNIa progenitors (see Fig. 1 of Mannucci 
2008). The integrated yields are given by: 

YUsnia = A bi „Pi(snia) IMF(M B )dM B , (6) 

where Mb is the mass of the binary system. y, jS ,„ fl values for ten 
/-elements and two different M up values are presented in Tables 
4 and 5. 

5) The evolution of all the models is followed from an 
initial time, fo=0, to a final time, ?y = 13.5Gyr, with a step 
of Af=0.01Gyr. The initial conditions in all the models are 
M gas (to)=0 and M. s , fl ,.(t )=0. 

6) For the outflow rates, we use two types of prescriptions: 

6.1) Well-mixed winds. We assume that the ejecta of 
MS, LIMS, and SNIa are well-mixed with the ISM of the galaxy 
before part of the ISM is expelled to the intergalactic medium 
(IGM). In this work, we consider that well mixed winds depend 
on the SFR and consequently the outflow rate is given by, 

w(t) = v«A(f), (7) 

where v is a free parameter in our models, adjusted to reproduce 
the observed M gas of the galaxy. 

6.2) Selective winds. We assume that during the time 
interval when a selective wind lasts, the ejecta produced by MS 
during that interval is expelled from the galaxy to the IGM with- 
out mixing with the ISM, therefore these ejecta do not contribute 
to the chemical enrichment of NGC 6822. The total mass lost 
due to the selective winds is small relative to the gaseous mass 
of the galaxy. 

4. Models for abundances derived from collisionally 
excited lines 

4.1. Observational constraints 

NGC 6822 has been well studied along the years by differ- 
ent authors, therefore good observational constraints can be 
found in the literature. For instance, CCPO derived its total dy- 
namic, total baryonic, and gaseous masses. They amount to 
M T = (2.0+0.5)xl0 10 Mq, M bar = (4.3.+0.2)xl0 8 M G and 
M,; flJ =(1.98±0.2)xl0 8 Mq, respectively. 



The chemical abundances of several H n regions, determined 
by diverse authors from collisionally excited lines (not corrected 
for dust depletion) were compiled by HPCG09, with the ex- 
ception of Cl/H that comes from Peimbert et al. (2005). The 
average O/H abundance of those Hn regions is 12+log(0/H)= 
8.08 + 0.05, and HPCG09 found no evidence of chemical inho- 
mogeneities in this galaxy for the H n regions located in a central 
area of about 3 kpc in diameter. In this work we will include a 
correction due to depletion in dust grains, by adding 0.08 dex to 
the O gaseous abundances of H n regions, correction suggested 
by Peimbert et al. (2005) as adequate for NGC 6822. From this 
we obtain 12+log(0/H)= 8.16 ± 0.05. On the other hand, from 
A-type supergiant stars, Venn et al. (2001) derived the chemical 
abundance of iron to be 12 + log(Fe/H)=7.01±0.20. These data 
will be used as observational constraints for the present compo- 
nent of the ISM. 

In addition, HPCG09 obtained the chemical abundances of 
1 1 PNe. They found that some of them are relatively young and 
their average O/H ratio (12+log(0/H)= 8.11 ± 0.10) is similar 
to the one of Hn regions. We tag these PNe as 'young PNe'. 
Besides, there is a sample of PNe with much lower O, showing 
an average of 12 + log(0/H) = 7.72+0. 10. We tag these PNe as 
'old PNe'. The chemical abundances of both groups will be used 
in the chemical evolution models, as additional constraints cor- 
responding to the past and recent-past component of the ISM. 
No correction for dust is considered for PNe. 

At the bottom of Table [Qwe present the O/H, N/O, Ne/O, 
S/O, Cl/O, and Ar/O abundance ratio averages for the best de- 
termined H n regions (the correction for dust is included), young 
PNe and old PNe. There is no observational value for C/H ob- 
tained from CELs, so there is no observational constraint for 
C/O. Following Allen et al. (1998), we will assume the follow- 
ing ages for the two PN groups: between 0.1 and 1 Gyr for the 
young PN population and between 3 and 9 Gyr for the old one. 
The present age of the galaxy is assumed to be 13.5 Gyr. Also 
in this table we present the Fe/H value by Venn et al. (2001) 
from A-type supergiants. Since most of the Fe in Hn regions 
and PNe is trapped in dust grains (Peimbert & Peimbert 2010; 
Shields 1978; Delgado-Inglada et al. 2009) we do not present 
the gaseous Fe/H abundance. In this table we are not showing 
either He, C and N abundances for PNe because PN progenitors 
modify strongly these element abundances during their evolu- 
tion, thus He, C and N abundances in the PN envelopes are not 
representative of those values in the ISM at the formation time 
of PN progenitors and cannot be considered as constraints for 
a chemical evolution model. Moreover, CI abundances of PNe 
are not shown in Table 1 , because they were not determined by 
HPCG09. 

4.2. The computed models 

We computed four chemical evolution models using different 
upper-mass limits for the IMF and different galactic winds pre- 
scriptions in order to reproduce the observed M gas and the O/H 
values (corrected for dust depletion) obtained from CELs in H n 
regions. We labeled them MIC to M4C. In Table Q] we present 
the characteristics and time evolution of these four models. In 
column 1 we list the model name; in column 2, the upper-mass 
limit for the IMF used in the model; in column 3, the correspond- 
ing wind prescription: 'W' for well-mixed and 'S' for selective 
wind; in column 4, the gaseous mass predicted by the model at 
the present time (i.e., 13.5 Gyr); in column 5, we list the time at 
which the chemical abundances are obtained, they correspond to 
the times adopted for H n regions, young PNe, and old PNe; and 
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in columns from 6 to 13, the predicted abundance ratios at the 
different times. As we mentioned in §4. 1, in the three rows at the 
bottom of this table, we include the observed values for H n re- 
gions, young PNe and old PNe, which are used as observational 
constraints. 

For model MIC we adopted an M up of 60 M Q for the IMF 
(equal to the value adopted by CCP06 in their best model), and 
we did not include any outflow. From Table Q]it is evident that 
MIC grossly fails to fit M gas as it predicts a gas mass ~4.5 times 
higher than the observed value. 

Thus, to reduce the enormous difference between the M gas 
predicted by model MIC and the observed value, we have com- 
puted model M2C with a well-mixed outflow during the early 
history of the galaxy. This was originally proposed by CCP06 
who computed the evolution of gaseous thermal energy in the 
galaxy provided by the SFR, under the assumption of no gas 
outflow (see their Fig. 7). They demonstrated that the accumu- 
lated thermal energy is much larger than the present binding en- 
ergy, therefore much of that energy should have been thrown 
away through galactic winds. CCP06 gave several arguments to 
explain why the outflows, needed to eliminate this excess in ther- 
mal energy, likely occurred in the first few billion years of evo- 
lution (t < 5 Gyr) when the galaxy was less bounded (note that 
the binding energy of a galaxy is also a function of time and in 
the hierarchical cosmology, the potential wells of galaxies for a 
given present mass, were shallower in the past). The gas outflow 
should have stopped when the binding and thermal energy came 
to an equilibrium. 

In our model, the well-mixed galactic wind is turned-on at 
t=1.2 Gyr, when the star formation starts and it lasts 5.3 Gyr, the 
time necessary for the ejected mass to be the maximum possible, 
leaving the minimum amount of gas required to maintain the 
SFR. During the wind phase, a value v = 6.67 in the outflow 
rate equation (Eq. [7]) was required. The evolution of the total 
baryonic mass and the gaseous mass in the galaxy, as a function 
of time, is shown in Fig. [2] M llp These to quantities are directly 
related with the infall, SFR and galactic wind. The dashed line 
is the time evolution of the M gas (Ec. 1) and the solid line in the 
time evolution of the total baryonic mass (dMt, ar /dt = f(t) - 
w(0). 
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t (Gyr) 

Fig. 2. Evolution of the total baryonic mass evolution (solid line) 
and gaseous mass (dashed line) used in our models M2C, M3C, 
M4C, and MIR. A well-mixed galactic wind occurs between 1.2 
and 6.5 Gyr. Observational constraints are as in CCP06. 



The present-time O/H value predicted by model M2C is 4 
cr higher than the observed values in Hn regions. Similarly the 
ISM O/H ratio found at f=12.9 Gyr is 2 cr higher than the ob- 
served values of the young PN population, and for the old PN 
population, the predicted ISM O/H value at ?=7.5 Gyr is 2.6 <x 
higher than observed (see Table Q]) 

Therefore, to better reproduce the O/H abundance ratios we 
constructed model M3C, by adding selective winds to the as- 
sumptions of model M2C. A selective wind, where only the ma- 
terial produced by MS is expelled to the IGM, might have been 
caused by the presence of a notable increase in the star formation 
rate in the last 600 Myr, as reported by CCP06 and Gouliermis et 
al. (2010, and references therein). Model M3C required a selec- 
tive wind between 12.5 and 13.5 Gyr in order to reproduce the 
current O/H value given by H n regions. Therefore this model is 
identical to model M2C for t<12 Gyr, but reproduces well the 
observed O/H ratio and it also reproduces well the O/H values 
for both PN populations. For the other elements which are pro- 
duced by MS (Ne, Ar, CI and S) the predicted Xi/O ratios are 
similar in both models, unlike those elements partially produced 
by LIMS (C and N), for which the abundance ratios relative to 
O are slightly larger than in M2C due to the loss of metals pro- 
duced by MS in the selective wind phase (see TableQ]). Although 
model M3C succeeds in fitting the observed O/H abundances, 
we do not consider it a satisfactory model because the physical 
mechanism to produce a 1 Gyr selective wind, where all chem- 
ical elements ejected by all MS are expelled to the IGM before 
they can mix with the ISM, is unlikely. 

Then, to fit the O/H ratios observed in H n regions, we pro- 
posed a lower value for M up , keeping the other assumptions of 
model M2C (a similar model was proposed by HPCG09). For 
model M4C we used an M up = 40 M , including a well-mixed 
wind, starting at f=1.2 Gyr and lasting 5.3 Gyr. Figure [3] shows 
the evolution predicted by model M4C for O, C, N, Ne, S, CI, Ar 
and Fe, relative to H, as a function of time. The thick lines (in the 
eight panels) represent the model results, the symbols represent 
the average values observed for H n regions and the two PN pop- 
ulations as explained in the figure caption. We consider model 
M4C as our best model to reproduce the chemical abundances 
derived from CELs, therefore we will discuss its characteristics 
in the following subsection. 

4.3. Our best model for abundances from CELs 

Model M4C was tailored to fit the observed O/H abundance ratio 
by reducing M up to 40 M . We find that it also reproduces well 
the present abundances of the elements heavier than O. This M up 
may be considered a low value for an IMF, however it seems a 
plausible value since some authors (e.g., Goodwin & Pagel 2005 
and references therein) have argued that dwarf galaxies could 
have smaller M up values than that of the solar vicinity. 

From the time evolution shown in Fig. [3] it is found that for 
all the elements, the temporal enrichment shows two minima, 
at 4 and 1 1 Gyr, due to the dilution caused by the large infall 
and to the decrease of the SFR (see Fig. [TJ, respectively. We 
find that the predicted abundances of all the elements at t ~ 7 
Gyr agrees with the average abundances of old PNe, whose pro- 
genitors were formed at t = 7.5 + 3.0 Gyr (6.0 +3.0 Gyr old). 
This epoch coincides with the beginning of the second fast en- 
richment. Furthermore, the recent abundances predicted for all 
the elements matches within uncertainties (~0.1 dex) the aver- 
age abundances of young PNe, whose progenitors were formed 
at t = 12.9 ± 0.5 Gyr (0.6 +0.5 Gyr old). In addition, Fig.[3]shows 
that model M4C matches well the H n region average data of all 
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elements within ~1 .5 a (which corresponds to the uncertainties 
in abundance determinations), except in the case of N where it 
fails grossly; this problem is discussed in more detail in the next 
paragraph. Finally, to fit the observed values of Fe/H, we needed 
to assume a fraction A^,-,, = 0.01. 

We want to remark that the model fails considerably in re- 
producing the presently observed N/H abundance ratio. This 
is found in several works on chemical evolution models (e.g., 
Chiappini et al. 2003; Carigi et al. 2005; Molla et al. 2006; 
Romano et al. 2010; Kobayashi et al. 2011). The prediction of 
our model is 7.4 <x higher than the value observed in Hn re- 
gions, therefore we consider that the N yields for LIMS adopted, 
in general, for the models have been overestimated. 

In order to confirm the results given in the previous para- 
graphs, in Fig. g] we present the C/O, N/O, Ne/O, S/O, Cl/O, 
Ar/O and Fe/O abundance ratios, as a function of 12+log(0/H), 
predicted by model M4C. Here we show the data for the best ob- 
served objects published by HPCG09. H n regions are presented 
as magenta filled circles, young-PNe are represented by green 
filled squares and old-PNe, by blue filled triangles. The open 
squares are the two type-I PNe, and the open triangle is PN6 (in 
HPCG09 nomenclature) catalogued as an old-PNe based mainly 
on its Ar/H abundance. We want to remark that in this figure we 
present the most complete set of observational constraints for 
NGC 6822 to be compared with our chemical evolution model. 

In this figure, regarding N/O, we find again that the model 
predicts to much N (probably due to problems with N yields). 
Unfortunately, we cannot comment about the C/H abundance be- 
cause there are no determinations of this element using CELs. 

In the case of Ne, all Ne/O ratios in PNe are above theH n 
values and above model predictions. This could be due to the 
ionization correction factor used to derive Ne (Ne/O = Ne ++ /0 ++ 
is commonly used) which is not adequate for low ionization H n 
regions or low density PNe (Peimbert et al. 1992, 1995b). Ne/O 
from the model agrees better (within uncertainties) with the val- 
ues for H ii regions which means that the Ne and O yields for 
massive stars up to 40 M Q , assumed in our model, are reliable. 

For S/O, the comparison is not good. In particular, the model 
predictions are 2 cr higher than the observed values in young 
PNe and Hn regions. The uncertainties in S abundance deter- 
minations are very large, in particular for PNe, therefore it is 
difficult to draw any conclusion, regarding S yields. 

For Cl/O our model is in good agreement with the present 
value of the ISM within the uncertainties. Similarly the evolution 
of Ar/O agrees, within the uncertainties, with the observational 
values. 

We reproduce the present Fe/H value obtained from A-type 
supergiant stars by assuming that only 1% of all the stars be- 
tween 3M and 15M become SNIa. 



5. Models for abundances derived from 
recombination lines 

It is well known that the C ++ and ++ abundances determined 
from the faint C n and O n recombination lines in photoionized 
nebulae are systematically larger than abundances derived from 
collisionally excited lines of the same ion (e.g., Rola & Stasinska 
1994; Peimbert et al. 1995; Liu et al. 2004; Esteban et al. 2004; 
Tsamis et al. 2004; Wesson et al. 2005; Peimbert & Peimbert 
2011; Simon-Diaz & Stasinska 2011; Tsamis et al. 2011; Carigi 
& Peimbert 2011). Other ionic abundances such as + /H + and 
N + /H + , derived from recombination lines show the same behav- 
ior. This problem is known as the 'Abundance Discrepancy" and 



is generally parametrized by the abundance discrepancy factor 
(ADF) defined as: 

ADF(X i+ ) = (X i+ /H + ) RL /(X i+ /H + )cEL. (8) 

ADFs have commonly values of about 2. Then, trying to discrim- 
inate between the abundance ratios obtained from both, CELs 
and RLs, we have computed a model to reproduce the O/H abun- 
dance calculated from RLs and we compared it with our best 
model based on collisionally excited lines (model M4C). 

5.1. Observational constraints 

We used the abundances reported by Peimbert et al. (2005), 
from RLs, for the Hn region HV. They obtained a value 12 
+ log(O/H)=8.42±0.06, where this determination includes the 
correction for dust depletion. Peimbert et al. (2005) considered 
dust corrections of 0.08 dex for O and 0.10 for C, these correc- 
tions have been adopted for H n regions throughout this paper. 
Abundance values for other elements are taken from Table 9 of 
Peimbert et al. (2005). 

On the other hand, there is no abundance determinations 
based on RLs for PNe in NGC 6822 because these faint lines 
were not detected in these objects. Thus we have not direct ob- 
servational constraints based on RLs for the past component of 
the ISM. Therefore we have used the determinations of chemi- 
cal abundances derived from CELs as reported by HPCG09 and 
we have modified them by considering a certain constant value 
for the ADF. Due to this procedure, the RLs abundances for PNe 
should be considered only indicative. 

Liu et al. (2006) report that for most galactic PNe, ADFs 
are in the range of 1.6 to 3. Likewise Peimbert et al. (1995a) 
determined ADFs for many galactic PNe and computed similar 
ADF values. Therefore we have applied an ADF=2 to the PN 
abundance determinations from CELs to derive a value for re- 
combination lines. Note that the ratio of any heavy element to 
O in PNe, is not affected by the ADF correction because the 
abundances of both elements have been increased by the same 
amount. 

In the bottom part of Table |2]we present the chemical abun- 
dances for H ii regions obtained via recombination lines plus cor- 
rection due to dust depletion. In addition we present the chemical 
abundances for the young and old PN populations calculated af- 
ter using the ADF mentioned above for the abundance values by 
HPCG09. 

5.2. The model 

We computed a single chemical evolution model using similar 
assumptions to those described for model M4C. The model is 
labeled MIR and it was tailored to reproduce the O/H abundance 
obtained for the H n region HV, from recombination lines. 

The same as in model M4C, model MIR required a well- 
mixed wind, starting at 1.2 Gyr and lasting 5.3 Gyr, to reproduce 
M gas . Now an M„ ; ,=80 M is needed to reproduce the current 
O/H value. In Table [2] (top part) we present the gaseous mass 
and the O/H, C/O, N/O, Ne/O, S/O, Cl/O, Ar/O, and Fe/H ratios 
predicted by model MIR for different times. 

In this model, the computed O/H value at t = 7.5 Gyr is in 
very good agreement with the adopted value for the old PN pop- 
ulation. For the young PN population, the oxygen abundance cal- 
culated, is also in very good agreement with the adopted value 
which is very similar to the one in Hn regions. To fit the ob- 
served Fe/H value we required to adopt, again, At,j„ = 0.01, 
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Fig. 3. Chemical evolution of elements in NGC 6822 as predicted by model M4C. To reproduce the observed M gas value and the 
different heavy element abundances obtained from collisionally excited lines for H n regions, a well-mixed galactic wind lasting 
from 1.2 to 6.5 Gyr and an M up -AQ M are assumed in this model. The O/H value of Hn regions has been corrected for dust 
depletion, as explained in §4. 1 . The filled circles, filled squares, and filled triangles represent the average observational values for 
H ii regions, young, and old PN populations, respectively. The filled star in the Fe/H panel represents the value derived by Venn et 
al. (2001) for A-type supergiant stars. Observed C/H value is not shown because there is no observed value from CELs. 



Fig. shows the time evolution of the element abundances 
(O, C, N, Ne, S, CI, Ar, and Fe) relative to H as predicted by 
model MIR. The thick lines in every panel represent our model 
results and the symbols are equal to those in Fig. 3. In this model 
the computed C/H value agrees well with the observational value 
obtained from RLs. The predicted N/H ratio is slightly closer to 
the observational constraint than in the model M4C, but it is still 
3.6 cr higher than the observed value, again indicating that the N 
yields for LIMS have to be revised. The predicted evolution of 
Ne/H shows relatively good agreement with observed values of 
old PNe but not so good for young PNe and Hn regions. In the 
case of S/H, the predicted evolution is in a good agreement with 
the PNe observed values but not with observations of Hn re- 
gions. The predictions for Cl/H and Ar/H are significantly lower 
than the observed values, in all the cases. 

The described behaviors are clearly noticed in Fig.|6]where 
we present the predicted evolution of the C/O, N/O, Ne/O, S/O, 
Cl/O, Ar/O and Fe/O abundance ratios as a function of 12+log 
O/H and the individual values of the best observed objects (see 
§4.3). In this figure, many of the Ne/O observed values (in par- 
ticular those of young PNe and H n regions) lie above the predic- 
tions, the S/O for Hn regions is well predicted within l<x level. 



The predicted Cl/O and Ar/O values are ~ 2 cr lower than the ob- 
served values. This latter problem is due to the flat extrapolation 
assumed for the SN yields for m > 40 M (see §3). According 
to Woosley & Weaver (1995), CI and Ar yields for m = 40 M 
are lower than for m = 35 M and therefore very massive stars 
would contribute less than massive stars to the ISM enrichment 
of those elements. 

We consider that to make a better comparison between mod- 
els and observations SN yields for m > 40 M are necessary. 
Recall that the yields in this paper for m > 40 M were extrap- 
olated from the m = 40 M yields and this extrapolation seems 
not adequate. 

Finally, to fit the observed log [Fe/H] =-0.49 value we 
adopted Ai,,„=0.01, the same value used for CELs models. 

5.3. The stellar yields problem 

One of the most important elements in the chemical evolution 
models are the stellar yields. In the literature there are several 
works about stellar evolution models that present stellar yields 
computed from different codes. These stellar models include dif- 
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Fig. 4. Element abundance ratios relative to O as predicted for model M4C compared to observed values of PNe and H n regions as 
presented in HPCG09. The O/H abundance of Hn regions has been corrected for dust depletion, as explained in §4.1. Hn regions 
are represented by magenta filled circles, young-PNe by green filled squares and old-PNe, by blue filled triangles. The open squares 
are the two type-I PNe, and the open triangle is PN6 (see text or HPCG09). The magenta filled star in the Fe/O panel represents the 
value derived by Venn et al. (2001) for A-type supergiants. Observed C/O value is not shown because there is no observational C/H 
value from CELs. 



ferent physical ingredients and it is difficult to compare them 
(e.g., Romano et al. 2010 and references therein). In addition, 
there is no complete grid of stellar yields available for the entire 
range of stellar masses. 

Our main problem to reproduce RLs observations for some 
chemical species is the lack of a homogeneous and complete 
set of stellar yields. First of all, as we explained in detail in §3, 
for LIMS we use the Karakas & Lattanzio (2007) yields, which 
cover from 1 to 6 M and for MS, those from the Geneva group 
(from 9 to 80 M Q ) for the pre-SN stage plus those from Woosley 
& Weaver (1995, hereafter WW) (for 9 to 40 M Q ) for the SN 
stage, together with the special treatment proposed by Carigi & 
Hernandez (2008). In order to reproduce the abundances from 
RLs, we need to extrapolate the WW yields from 40 M to 80 
M . We can see in Figs. 3 and 5 that the predictions for the evo- 
lution of O/H and Ne/H are in good agreement with the obser- 
vational constraints obtained from CELs (model M4C) and from 
RLs (model MIR) for the present (Hn regions) and past (PNe) 
components of the ISM. However for the other chemical ele- 
ments, C/H, N/H, S/H, Cl/H, Ar/H, and Fe/H the predictions do 
not change for the two models, even when for model MIR M up 
is 80 M Q . 



Moreover, Cl/H and Ar/H values predicted by model M4C 
are in good agreement with observations, but this is not the case 
for model MIR. We blame this to the lack of SN yields for 
masses > 40 M Q and the extrapolation made by us up to 80 M Q . 
We have to keep in mind that M„ ; , of the model M4C is equal 
to 40 M and no extrapolation for the yields was needed, while 
M up of MIR is equal to 80 M and we have to extrapolate the 
yields. 

Once we have these results, we did an exercise changing 
only the stellar yields of massive stars. We took the SN yields by 
Kobayashi et al. (2006) and the pre-SN yields by the Geneva's 
group. We kept the rest of assumptions described in §3 and 
constructed two models, M4Ckoba (solid line) and MIRkoba 
(dashed line) presented in Fig Q 

We plot the evolution of the O/H, Ne/H, S/H, Cl/H, and Ar/H 
in order to compare them with the observations. The predictions 
for the evolution of Ne/H and S/H abundances are in very good 
agreement with the observational constraints obtained from RLs 
only (dashed line and open symbols) but for CELs (solid line and 
filled symbols) the model predictions are ~ 3<x higher. Based on 
O/H, Ne/H and S/H, Kobayashi et al. (2006) yields produce a 
good model for observed RLs abundances but not for observed 
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Fig. 5. Chemical evolution of elements in NGC 6822 as predicted by model MIR. To reproduce the observed M gas value and the 
different heavy element abundance ratios obtained from RLs for Hn regions, a well-mixed galactic wind lasting from 1.2 to 6.5 
Gyr and an M, (; ,=80 M are assumed in this model. The observed O/H and C/H values have been corrected for dust depletion, as 
explained in §4. 1 . The filled circles, filled squares, and filled triangles represent the average observational values for H n regions, 
young, and old PN populations, respectively. The filled star in the Fe/H panel represents the value derived by Venn et al. (2001) for 
A-type supergiant stars. 



CELs abundances. If we want to force the M4Ckoba model to 
reproduce abundances from CELs, we must diminish even more 
M up in the IMF (M4C has an M„ p =4OM ). The Cl/H and Ar/H 
predictions of both models are lower than the observational con- 
straints, by ~ 0.75 and ~ 0.25 dex, respectively. In order to re- 
produce the Cl/H and Ar/H observed abundances in H n regions, 
the CI and Ar yields should be 5.6 and 1.8 times higher than 
those computed by Kobayashi et al. (2006). This is an impor- 
tant result, because previous works that test stellar yields (e.g., 
Timmes et al 1995; Romano et al. 2010; Kobayashi et al. 201 1) 
cannot check the CI and Ar yields due to the lack of abundance 
determinations of these elements in stars of the solar vicinity. 



Based on this exercise, we prefer to keep the SN yields by 
WW for the model reproducing abundances from CELs (M up = 
40 M ). Nevertheless, we prefer the SN yields by Kobayashi et 
al. (2006) for the model reproducing RLs abundances (M llp = 80 
M ) with CI and Ar yields increased by factors of ~ 6 and ~ 2, 
respectively. 



6. The AY/AO ratio and the IMF 

Our chemical models predict the evolution of the He and O en- 
richment, both by mass, known as AY and AO, as a function of 
time. For convenience, and to compare with observations the AY 
/AO ratio is often used, where AO is equal to O and AY is equal 
to Y - Y p . Here Y p is the primordial He abundance and amounts 
to 0.248, the value derived from metal poor irregular galaxies 
and WMAP (Peimbert et al. 2007a; Dunkley et al. 2009). It is 
well known that AY versus AO mainly depends on the adopted 
yields and the IMF. 

Fig. [8] shows the AY vs. AO behavior for our models M4C 
(solid line, this model reproduces the O/H abundances given by 
CELs) and MIR (dashed line, reproducing the O/H abundances 
given by RLs ). Both lines show a similar behavior (same form, 
different slope) which depends mainly on the history of the mass 
accretion, the variation of the SFR with time, and on the time de- 
lays assigned to the LIMS for their contribution to the chemical 
enrichment of He. 

In Fig. [8] the initial slope corresponds to the first epochs of 
the galaxy, where only infall and IRA from MS are the processes 
controlling the He and O abundances. The sudden discontinuity 
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Fig. 6. Element abundance ratios, relative to O, as predicted for model MIR compared to observed values from recombination 
lines for the Hn region HV as presented in Peimbert et al. (2005) and the adopted values for PNe, as explained in §5.1. O and C 
abundances in HV have been corrected for dust depletion, as explained in §4. 1 . H n region HV is presented in magenta filled circl, 
young-PNe are represented by green filled squares and old-PNe, by blue filled triangles. The open squares are the two type-I PNe 
and the open triangle is PN6 (see HPCG09). The magenta filled star, in the Fe/O panel, represents the value derived by Venn et al. 
(2001) for A-type supergiant stars. 



with a large increase in He, occurs when LIMS start their contri- 
bution to the He enrichment of the ISM; these stars mainly pro- 
duce He but not O. From this point onwards both curves evolve 
in a continuous increase of He and O, with an almost constant 
slope, perturbed only by a small loop which occurs when the 
SFR changes suddenly, in the epoch from 6.3 to 10.5 Gyr (see 
Fig- CD bottom). At even later times both curves increase with an 
apparently constant slope. 

Interestingly, although the behavior is similar, both curves 
showing the same perturbations due to the behavior of the SFR, 
the slopes are very different. As expected, the slope is larger for 
M4C, which has M up of 40 M Q , and thus in this model less oxy- 
gen is produced. The slope derived for model M4C, considering 
only the zone from AO > 0.5 x 10~ 3 , is about 7.2, while for 
model MIR it is about 3.7. 

In Fig.[8]we also include two observational points. The filled 
square represents the (AY, AO) values derived from the He i RLs 
by Peimbert et al. (2005) and the [Om] CEL (HPCG09) un- 
der the assumption of a constant temperature given by the [O m] 
CELs. The open square represents the values derived from abun- 
dances obtained through He i and O n RLs, considering the pres- 
ence of spatial temperature variations. It is observed that, in both 



cases, AO is well fitted because the models were tailored to fit 
O/H. On the other hand, the observed AY shows huge uncertain- 
ties, because it is very difficult to derive a precise value for He 
abundances due to: i) the large statistical errors due to the faint- 
ness of the lines (represented in the Fig. [8] by the error bars), 
and ii) the possible presence of neutral He in the observed H n 
regions, that would increase the derived AY value for both deter- 
minations. Both Y values are derived from RLs, one under the 
assumption of constant temperature and the other under the as- 
sumption of the presence of temperature variations, the second 
one gives a Y value smaller than the primordial value probably 
indicating the presence of neutral helium inside the observed H n 
region. 

We can compare also the slopes predicted by our models 
with observational data in the literature. For instance, Izotov et 
al. (2007) have derived chemical abundances for a large number 
of Hn galaxies of different metallicities. From these data they 
derived a slope AY/AO = 3.6+0.7 or 3.3 +0.6, depending on the 
He i emissivities used, these numbers have not been corrected by 
the fraction of O atoms trapped in dust grains. By considering 
that such dust fraction amounts to 0.08 dex (Peimbert et al. 2005) 
we obtain from Izotov et al. (2007) data that AY/AO = 3.0+0.6 
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Fig. 7. We present the evolution for O/H, Ne/H, S/H Cl/H and 
Ar/H for two models, M4Ckoba (solid line) and MIRkoba 
(dashed line), based on Kobayashi et al. (2006) yields. The filled 
circles, filled squares, and filled triangles represent the average 
observational values for H n regions, young, and old PN popula- 
tions, respectively. 
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Fig. 8. Evolution of AY vs. AO for our two best models, M4C 
(solid line) and MIR (dashed line). The black square represents 
observed abundances under the assumption of constant tempera- 
ture, and the open square represents abundances from RLs under 
the assumption of temperature variations (Peimbert et al. 2005), 
note that the observed AY values might be lower limits due to the 
possible presence of neutral helium inside the H n region. 



or 2.75 +0.5. On the other hand, Peimbert et al. (2007b), from 
theoretical and observational results, obtain AY/AO = 3.3±0.7. 
These slopes are in better agreement with model MIR, which 
was computed to fit O abundances from RL and for which an 
M up = 80 M was assumed. 

To summarize, model MIR produces a AY/ AO slope similar 
to the one derived for a set of other galaxies, while M4C does 
not. In the case of NGC 6822, however, a much better Y determi- 
nation is needed to discuss this problem further. 



7. Discussion and conclusions 

We have developed a code to calculate chemical evolution mod- 
els for galaxies. The code includes the instant recycling approx- 
imation for MS, and a time delay for the contribution of LIMS. 
Such a time delay has been obtained for each element following 
the prescription given in Franco & Carigi (2008). Our code also 
includes the SNIa contribution to the chemical evolution in the 
galaxy. The code can be used for different SFRs, IMFs, yields, 
infalls and outflows. It predicts the evolution of 27 elements. 

Using this code we calculated chemical evolution models for 
the dwarf irregular NGC 6822. The main characteristics of the 
models for this galaxy, are: 

- A time dependent infall, parametrized from the one obtained 
by CCP06 from cosmological considerations, with a primordial 
abundance given by X=0.752, F=0.248 and Z=0.000 (Fig.Q}. 

- The star formation history derived by CCP06, which is based 
on the photometric properties of the galaxy (FigQ}. 

- The IMF by Kroupa et al. (1993) with different M up values 
depending on the model. 

- A set of metallicity dependent yields for 27 elements including 
H, He, C, N, O, Ne, S, CI, Ar, Fe and others. 

- 1% of the stars with masses between 3 and 15 M gives up 
binary systems and every one of these systems becomes a SNIa. 
The contribution of SNIa enrich the ISM with a time delay of 
1 Gyr after the formation of the SNIa precursor. 

- Two types of time dependent outflows can be used: well mixed 
and selective. 

Our models for NGC 6822 were constructed to fit the follow- 
ing observational constraints: a) M gas , b) O/H chemical abun- 
dances of Hn regions and the Fe/H ratio of A-type supergiant 
stars representing the current abundances of the ISM, and c) 
chemical abundances of young (0.6+0.5 Gyr) and old (6±3 Gyr) 
planetary nebulae representing the abundances of the ISM at the 
time they were formed. The H n abundances have been corrected 
for the fraction of C and O in dust grains. Two sets of chemical 
abundances were used: those derived from CELs (HPCG09) and 
those derived from RLs (Peimbert et al. 2005), they differ by 
an abundance discrepancy factor, ADF, of about 0.26 dex. As 
the RLs were not observed in PNe, we estimated their RL abun- 
dances based on the CEL abundances, and adopting an ADF of 
2. Thus PNe abundances from RLs are only indicative. 

The baryonic total mass captured by the system is given 
by the gas inflow estimated from cosmological considerations, 
and the amount of gas transformed into stellar mass is obtained 
from photometric data. From these two considerations only, the 
remaining gas turns out to be 4.5 times higher than observed. 
Therefore to reproduce the observed M gas value the models re- 
quired a well mixed outflow lasting from t — 1.2 Gyr (which is 
the time when the star formation starts) to t — 6.5 Gyr (when the 
gas mass has diminished enough to fit the present gaseous mass 
but without interrupting the star formation, Fig.|2|. 

The SFR used in our models was derived from photometric 
data by CCP06, and it is in good agreement with the abundances 
obtained for old PNe (HPCG09). These two studies are inde- 
pendent and therefore old PNe data support the SFR based on 
photometric data. 

The two best models obtained to reproduce the present day 
ISM O/H value required: a) an M up = 40 M Q for model M4C, 
which is based on abundances derived from CELs, and b) an 
M up = 80 M for model MIR, which is based on abundances 
derived from RLs. 

Both models reproduce the O/H evolution of the past compo- 
nent of the ISM obtained from young and old PNe. This is par- 
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ticularly true for model M4C based on CELs. This result implies 
that the adopted SFR, derived from photometric data, describes 
well the shape of the O/H evolution and that the inferred well 
mixed wind is consistent with the absolute O/H value. 

The O/H ratio depends strongly on the use of RLs or CELs, 
but since the models are tailored to fit the observed values, it 
is not possible to decide which abundances are more adequate 
based on O/H alone. 

We have tested the time evolution of other elements (C, N, 
Ne, S, CI, Ar and Fe, relative to H). In general, the predictions 
of model M4C agree with the observations of H n regions and 
old and young PNe, indicating that the assumptions used for this 
model are correct and that the SN yields used in the 9 to 40 M 
range, derived from WW, are adequate. 

About 50% of C is produced by LIMS and about 50% by MS 
in the pre-SN phase. The C/O ratio is well fitted by model MIR, 
result that implies that the C and O yields are adequate. 

On the other hand the time evolution of heavy elements other 
than O, predicted by model MIR, are not always in agreement 
with the observations. Given that the predicted abundance ratios 
depend strongly on the yields, the adopted yields can then be 
tested. The C/O, S/O and Fe/O ratios predicted are in fair agree- 
ment with the observations. But for Ne/O, Cl/O and Ar/O, the 
predicted values are smaller than observed, probably indicating 
that the extrapolation of the SN yields from 40 to 80 M made 
in this paper is not good enough. A different set of yields for this 
mass range, provided by Kobayashi et al. (2006) was tested, find- 
ing a good agreement with observational values in both models 
(M4Ckoba and MIRkoba, Fig[7]) for Ne/H and S/H, but Cl/H and 
Ar/H are underestimated by factors of ~6 and ~2 respectively. It 
is worth to mention that this is the first time that observational 
constraints for these elements are presented and compared with 
the models in order to test the stellar yields. 

Model MIRkoba is in very good agreement with the ob- 
served value of O/H obtained from RLs, but in order to repro- 
duce the observed value obtained from CELs (model M4Ckoba) 
it is necessary to lower M up below 40 M Q . We require a new set 
of yields for the 40 to 80 M mass range to make a better test of 
the model MIR against observations. 

The predicted N/H values in both models are always con- 
siderably higher than observed and probably indicate that the N 
yields need improvement. About 80% of N is produced by LIMS 
and about 20% by MS in the pre-SN stage. Therefore lower N 
yields for LIMS would provide a better fit to the observed N/O 
values. 

At this point, we are not able to conclude which set of abun- 
dances (based on CELs or RLs) represents better the real abun- 
dances in the ISM. Model M4C , that assumes M up = 40 M , re- 
produces well the O/H, Ne/H, S/H, Ar/H, Cl/H, and Fe/H obser- 
vational constraints. While model MIR, that assumes M up = 80 
M , reproduces well the O/H, C/H, Ne/H, S/H, and Fe/H obser- 
vational constraints, but fails to reproduce Cl/H and Ar/H prob- 
ably due to inadequate SN yields for stars with masses above 40 
M . 

Since AY does not depend strongly on the adopted temper- 
ature and AO does, AY/ AO is a very good way to discriminate 
between MIR and M4C models. Unfortunately the helium abun- 
dance determination for NGC 6822 is poor due to the errors in 
the He I line intensities and the correction due to the neutral 
He inside the Hn regions. We have used the AY/ AO determina- 
tions derived by other groups, based on other galaxies, and we 
find that the MIR model agrees better with the observations of 
the other galaxies than the M4C model. A much more precise 



He abundance determination of the H n regions in NGC 6822 is 
needed to advance in this topic. 
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Table 1. Models that reproduce abundances calculated via collisional excited lines. The 'W and 'S' are for well-mixed and selective wind prescription, respectively. 
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a Abundances in units of 12 + log(0/H). 
* CCP06. 

c Adding 0.08 dex to O gaseous value due to depletion in dust grains. 
d Peimbert, Peimbert & Ruiz (2005). 

e Fe/H from A-stars by Venn et al. (2001) and CEL O/H from Hn regions. 



Table 2. Model that reproduces abundances calculated via recombination lines plus correction for dust depletion. In column 3, 'W is for well mixed wind prescription. 
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e Fe/H from A-stars by Venn et al. (2001) and RL O/H from Hn regions. 
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8. Appendix A 

In this section we present the returned masses, time delays, and 
integrated yields for the chemical elements considered in this 
work. 

In Table 3 we show the time delays of LIMS for the returned 
mass and for the integrated yields of eight chemical elements. 
These data are independent of the M up value adopted in the IMF, 
but they are dependent on the initial stellar metallicity Z,. 

The data shown in Tables 4 and 5 were obtained by assuming 
M up = 40 M and M up = 80 M , respectively. 
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Table 3. Time delay (Gyr) of LIMS. 
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Table 4. Returned mass and Integrated Yields for M up = 40 M . 
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" SNIa yields are independent of the initial metallicity. 



Table 5. Returned mass and Integrated Yields assuming M up = 80 M . 
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1.445x10" 


-5 


3.470x10-'° 


4.278xl0- 5 


6.043x10" 


-/ 


2.516X10- 3 


4.009x10" 


-8 


3.769x10- 


12 


1.834xl0- 4 



" SNIa yields are independent of the initial metallicity. 
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